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Abstract 



We determine the magnitude, phase, and s-dependence of p-uj "mixing" in 
the pion form factor in the time-hke region through fits to e~^e~ TT~^7r~ 
data. The associated systematic errors in these quantities, arising from the 
functional form used to fit the p resonance, are small. The systematic errors 
in the p mass and width, however, are larger than previously estimated. 
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I. INTRODUCTION 



The pion form factor Ft^{s) in the time-hke region is extracted from cross section mea- 
surements of e'^e~ —>■ 7r+7r~, namely, 

a{s) = aUs)\FAs)\\ (1) 

where crem('S) is the cross section to produce a structureless tt+tt^ pair and s is the usual 
Mandelstam variable. The isovector p resonance dominates the cross section for y/s ranging 
from 600 to 900 MeV, though isospin violation allows the isoscalar uj resonance to contribute 
as well. We wish to extract the rho-omega "mixing" matrix element Ilp;^{s), which to leading 
order in isospin violation is given by 



FAs) = Fpis) 1 + ^ 7^;^.^ , , (2) 




from data. Note that Fp{s) is the pion form factor in the absence of isospin violation. Fp{s) 
is subject only to general theoretical constraints; it is our purpose to determine how its non- 
uniqueness impacts the extraction of Ilpuj- Moreover, the p mass and width are themselves 
sensitive to the choice of Fp{s) and we wish to determine the systematic error in these 

quantities as well. 

Maltman et al. have discussed the separation of 11^^(5) into two contributions: one from 
the direct coupling oi oj ^ 2tx and the other from mixing, oj ^ p ^ 2-n [^^H. Such 
a separation is model dependent, and we shall not pursue it further. Rather, we wish to 
determine the constraint e^e~ h^tt' data places on the sum of these contributions; we 
term Yipt^{s) the effective mixing matrix element. Iipu,{s) is usually assumed to be both real 
and approximately s-independent 0; we wish also to test these assumptions in the context 
of e^e~ — > TT+TT^ data. An explicit s-dependence in Iip^{s) emerges as a consequence of 
the inclusion of a non- resonant contribution to — 27r 0,0; thus, the s-dependence of the 
effective mixing amplitude in the resonance region may partially constrain the role of this 
contribution. 

Previous determinations of 'n.pu){s) have used either the empirical uj ^ 271 branching ra- 
tio P,p!0| or e+e^ tt+tt^ data We prefer the latter method for several reasons. 
Information on the phase and s-dependence of np^(s) is not accessible from the uj ^ 2t[ 
branching ratio. Moreover, the determined p-uo mixing matrix element then explicitly de- 
pends on the relatively poorly known p resonance parameters. 

Our work differs from earlier analyses in that we enumerate a variety of forms for Fp{s) 
which satisfy the known theoretical constraints. We extract the magnitude, phase, and s- 
dependence of the effective p-uj mixing matrix element np^(s) from e^e^ tc^tc^ data and 
study the systematic error in the above parameters resulting from the choice of Fp{s). We 
have also studied the systematic errors in the p parameter extraction and find them to be 
much larger than previously reported. 

To extract Ilp^{s) from e"'"e^ tt+tt^ data we naturally wish to use the best data set 
available. The most recent data for the pion form factor in the time-like region is due to 



Barkov et al. [12|. They report nip = 775.9 ±0.8 ±0.8 MeV and r„ = 150.5 ±1.6 ±2.5 MeV, 



where the errors respectively refer to the error arising from the statistical and systematic 



2 



uncertainties in the experimental data and to the systematic error resulting from the form 
chosen for the pion form factor The Barkov et al. value for the p mass contributes 

some 40% of the in the Particle Data Group's 1994 world average for this quantity [jT3[ . 
In the 1996 compilation, the Particle Data Group rule the Barkov et al. p mass "probably 
wrong" on statistical grounds and exclude the Barkov et al. values from the world averages 
for the p mass and width [l^Jl|] . The only determination of the p mass currently included in 
the Particle Data Group world average based, at least in part, on time-like pion form factor 



data is that of Heyn and Lang |[T^. We apply our analysis to both the data included in 



Barkov et al. and to the time-like region data included in Heyn and Lang JT3| , p!7| , to 

ascertain whether any systematic differences exist between the data sets. 

The p resonance is relatively broad, so that the reported resonance parameters are nu- 
merically sensitive to the convention under which the mass and width are defined. The 
appearance of a resonance is associated with a complex pole at s = Sp in the elastic scat- 
tering amplitude, and this complex pole can be used to define the resonance's mass and 
width |l^. The separation of Sp into a mass and width is not unique, however, for both 



the real and imaginary parts of Sp appear in any physical process [|T^. Thus, whether one 
defines Sp = m? — imT, so that m = wRe Sp, or = m — iT/2, so that m = Ke^/s^ 



is a matter of convention. We will present results for the p mass and width under both 
conventions. 



II. THE PION FORM FACTOR AND p - uj "MIXING" 

Only general theoretical constraints guide the construction of the pion form factor in the 
time-like region. Charge conservation requires the form factor to be unity at s = 0: 

FM = 1 • (3) 

Moreover, it should be an analytic function in the complex s plane, with a branch cut along 
the real axis beginning at the two-pion threshold, s = 4m^. Finally, time- reversal invariance 
and the unitarity of the S matrix requires that the phase of the form factor be that of Z = 1, 
1 = 1 vr-vr scattering [^]. This last emerges as vr-vr scattering in the relevant channel is very 
nearly elastic from threshold through s ~ (rriTr + ^«^)^ ||2l| , p!5[| . In this region of s, then, the 



form factor is related to the / = 1, / = 1 phase shift. Si, via 

FAs) = e''''F:{s) (4) 

so that 

^1 , , Im Ft,(s) 

The above is a special case of what is sometimes called the Fermi- Watson-Aidzu phase 
theorem [gg,^. 

In the resonance region the phase and analyticity constraints can be realized via the 
Breit-Wigner form 

mlil + e) 

hm F^s) = — — , (6) 
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where £ is a real constant. The complex pole Sp associated with the appearance of the 



resonance is given by Sp = — inipT p. We have adopted the y Re Sp convention for the p 

mass. Alternatively, one could have written ^Js^ = rfip — iV p/2^ so that m = Rey^, but 
this is merely a matter of convention. Note that the m and fh = Re^^ prescriptions are 
related via fil 



m 



m 



r. (7) 



The Breit-Wigner form Eq. (j^) only satisfies the phase constraint as s 



a more general form for the pion form factor, suitable for all s, is needed. Generalizations 



satisfying the enumerated constraints have been constructed by various authors [24,15,25 



we will follow the work of Gounaris and Sakurai and of Heyn and Lang [|T5| in what 
follows. We include p-u mixing as per Eq. (D, so that our enumerated constraints are 
brought to bear on the form of Fp{s) alone, for the violations of the above constraints due 
to isospin breaking are small. 

Gounaris and Sakurai consider a Fp{s) of the form /(0)//(s) with /(s) such that 

f (s) = (e/^s) cot 61- tie/ ^s), (8) 

where /(O) is real. Note that both the normalization and phase constraints are manifest 
in such a construction. The 6l phase shift is parametrized via a generalized effective-range 
formula of Chew-Mandelstam type p6|, with two free parameters, a' and b', 

(P/v^) cot 6l = a' + b'k^ + k^h{s) , (9) 

where k and h{s) are chosen to be 

k = {s/4 - mlfl\ h{s) = ^ log (4^) ; ^ > 4< 



i(ml - s/4)^/^ h{s) = — ^arccot ; < s < Ami . (10) 

nJs \ Ami — s / 



In this manner the required analytic structure is imposed as well. The parameters a' and h' 
can be replaced by functions of and T p by noting the resonance conditions cot 5\ {m"^) = 
and 6l'{mVj = l/lmpTp), so that the resulting form factor is of the form |24 



F<^^(s) = -jml + dmpTp) 

' ^ ^ s-ml- Tp{ml/kl)[k^{h - hp) -is- m^k^^h'^] + tmpTp{s) ' ^ ^ 



where 



1 \ 3 

k \ mn 



kp = k{ml), hp = h{ml), r,(.)=rJ-l (12) 



and d is fixed in terms of mp and m^. 
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d 



3ml 



log 



nip + 2kp 



2m^ 



+ 



nip mp-m'-^ 



Tckp 



Note that the Breit-Wigner form, Eq. 



is recovered as s 



(13) 



It turns out that the 



two-parameter fit of Eq. (|Tl]) does not suffice to fit the pion form factor data in the p 
resonance region, and a modification of F^^(s) to include a multiplicative real function of 
s, constrained by only the normalization condition, is phenomenologically necessary. 

It is critical to note that the adoption of an s-dependent width, r(s), predicated by the 



phase constraint of Eq. 
of the complex pole Sp 



), implies that s 



im 



pT p no longer determines the position 



Sp = m — irUpTp, where f{sp) = 0, then for s 



Rather, if one were to determine the mass and width from 

the phase shift would nothe n/2, and 
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m 



the resultant form factor near s 

to determine the resonance parameters by comparing f{s) to s — — inipT p at the s in the 



would not be of Breit-Wigner form. We prefer, then, 

2 



physical region where the real part of f{s) vanishes. The resonance condition cot 5 
thus satisfied as s rn? . Our procedure is consistent with that of Ref. pi . 



Heyn and Lang write Fp{s) 



as 



7r/2 is 



(14) 



where for < s<(m^ -|- m,r)^, -^red(s) is purely real and can be approximated by a cubic 
polynomial [|l^, 

Fred(s) = 1+/91S + /?2S'+/?3S^ (15) 

All phase information is contained in the Omnes function VL{s) p8|, which is chosen to be 



n{s) 



Sp- s 



as"^ + bs + c — {s — 4m^)5f(s)/4 



Modulo the {sp — s)/sp factor, VL{s) also has a f{0)/f{s) structure, where now 

/(s) = as'^ + bs + c - {s ~ 'iml)g{s)/4: 
and g{s) is determined by the one-pion-loop diagrams in the p self-energy [pU|-pT 



(16) 



(17) 



1 — 4m^ /s ,s > 4m^ 



g[s) = u log \- lu ,u = 

TT 1 — U 

= u arctan — ,u = J Am? / s — 1 , < s < Ami 

TT M * 

1 , U + 1 I z— 

= u log ,u = \ 1 — Am-L s , s < 

TT U — 1 * 



(18) 



^(0) = -2/7r . 

The p resonance is associated with a zero of f{s) in the complex plane for s > Am^, yet f{s) 
also vanishes on the negative real axis when s = Sp, by definition. This zero simulates the 
left-hand cut in the N/D construction of the amplitude and physically corresponds to 



a bound state [^. The {sp — s)/sp factor in ^l(s) removes this singularity while preserving 
the normalization constraint. The singularity structure of Q{s) occurs in Fqs{s) as well; 
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the zero there, though, is at such large negative s that the inclusion of a (sp — s)/sp factor 



would be of no phenomenological impact |2J]. Two parameters are needed to describe the 



left-hand cut. The third free parameter in Eq. ( p!?] ) functions as a CDD parameter; the N/D 
solution of the partial-wave amplitude dispersion relation does not uniquely follow from the 
information input on the left-hand cut p3|j2^ . 

Once a, b, and c are determined from a fit to data, the p resonance parameters can also 
be determined. As discussed earlier in the case of the Gounaris and Sakurai form factor P3 



Eq. (|TTD, we determine the resonance parameters by comparing f{s) to s — — iiripT p at 
the s in the physical region where the real part of f{s) vanishes. This has the effect of 
requiring that f{s) be of Breit-Wigner form as s — > m^. Consequently, is determined by 
Re{/(m^)} = and Vp is determined by Im{/(m^)} = —nipTp HTsj . 

The structure chosen for f{s) in Eq. (|1^) is formally consistent with the phase constraint 
resulting from unitarity and time-reversal invariance, yet the constraint may not be numer- 
ically well-satisfied, as the parameters are fit to the time-like pion form factor data. We 
can, however, require that the parameters in f{s) reproduce the empirical / = !,/ = ! vr-vr 
scattering length a\ to gauge whether the results we extract for the p parameters and Ilp^ 
are sensitive to this additional constraint. We define Cl-^ clS 

\ = lim cot 6l , (19) 

CLi s^4m| 



where empirically a\ = (0.038 ± 0.002) m^^ [34|. Note that l/a\ is equal to 2m,ro' in 



the Gounaris and Sakurai model, Eq. (^, and is equal to (16m^a + 4m^6 + c)m^ in the 
Heyn-Lang model, Eq. (|16]). 



The pion form factor we fit to data is 



1 + 



1 np^(s 



3s- ml + im^T^ 



(20) 



We work to leading order in isospin violation, and we adopt the SU(3) value of 1/3 for the 
ratio of the electromagnetic coupling of the u to that of the p. Note that -Fred(s) is given by 
Eq. (|T5|) and Q{s) by Eq. (|l6l), though we will also replace fl{s) by Fqs{s), Eq. ([Tl|) , and by 
an Omnes function modified to resemble Fqs{s), 

^^^^'^ - bs + c-is-Aml)gis)/A ' ^^^^ 

We thus test the sensitivity of the fit to the specific manner in which the phase constraint 
is realized. Rather than Eq. (0), Barkov et al. |jT^ adopt a pion form factor such that 



I + + api + Upti 

where a^, a pi, and are real constants, to be determined from a fit to data, and Fy^(s), 
with V a vector meson, is determined by Eqs. (0,|r^,0) with nip, V p replaced by my, Ey. 
Fy^(s) is complex in the region of s we consider, so that Eq. (|22|) would seem to violate the 
phase constraint of Eq. (Rl). We fit Eq. (B^) as well to the time-like pion form factor data. 
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and we will examine the above issue through direct comparison with the / = !,/ = ! phase 
shifts. 

If the inelastic contributions to vr-vr scattering in the s region of interest are generated by 



isospin violation only, the solution of the Muskhelishvili-Omnes integral equation ||2^ with 
inelastic unitarity suggests that the effective p-u mixing matrix element, Ilpij{s), is real [0. 
It is also thought to be weakly s-dependent 0. Thus, unless otherwise stated, we make the 
replacement Ilp^{s) = Ilpuj{m'^). We can, however, also test these assumptions by replacing 
Ilp^{s) with 

IlpUs) ^ fi-U^l) + ^Im nUml) (23) 

or with 

%u{s) = flp^iml) + {s~ ml)%^{ml) , (24) 

where it is natural to expand the s-dependence of the effective mixing matrix element about 
s = m^. In this manner we can test whether the value of Re{Ilp^{m'^)} is sensitive to the 
matrix element's possible s dependence or phase. 



III. RESULTS 



We optimize the fit parameters using MINUIT [^. The parameters can be highly 
correlated, so that we compute the correlation coefficients and final errors with a double 
precision FORTRAN code using standard techniques [^. We fit data from threshold. 



= 2m^ ~ 280 MeV, through ^/s ~ + = 923 MeV, as this is the region in which 



the empirical / = 1, vr = 1 phase shift is essentially elastic |2T[]. Our data set consists of the 
82 points in this energy range included in the analysis of Barkov et al. |]T^,|TB|. The following 
fits are based on Eq. pO| ) — the forms chosen for f2(s), Fred(s), and Iipui{s) in each fit are 
indicated in Table |. For fits B and C above, a fit of the type indicated by A' and A^ has 
been performed as well, so that B' and B\ e.g., indicate B-type fits in which 'n.p^{s) has been 
replaced by Eq. (^) and Eq. (^31) , respectively. Note that in fit B the parameter c is fixed 
so that the model reproduces the empirical scattering length, a\ = (0.038 ± 0.002) pi| , 
as per Eq. (plOl) . For definiteness, note that we use = 139.57 MeV, = 781.94 MeV, 
and = 8.43 MeV in all fits. The results are given in Tables ^ and |IT1[ The given 
parameter errors arise from the optimization of the fit to data, whose errors include both 
statistical and experimental systematic uncertainties. Specifically, the parameter errors are 



given by the square root of the diagonal elements of the inverse of the curvature matrix [36 



The parameters are correlated, so that the errors in the "Output" of Tables |I| and |T| are 



generated using the full error matrix. All the parametrized forms are able to fit the data 
exceedingly well — x^/dof ~ 1 in all cases. Fits A and B, along with the Heyn-Lang A-type 



fit [jT5[ to the 1978 world data |T^, are plotted as a function of the pion-pair invariant mass 



q, q = i/s, with the data [|l^] included in the compilation of Ref. |]T2| in Fig. [T|. The shapes 
of A and B are very similar in the resonance region but are visibly different at small q as 
fit B is constrained to reproduce the empirical a\ scattering length. The scattering length 
extracted from fit A is nearly a factor of 2 larger than the empirical value; later we will 
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explicitly compare the phase of the fits we generate with the measured / = 1, / = 1 vr-vr 
phase shift in the s range of interest. 

The effective p-ut mixing matrix element Ilpt^(m^) is remarkably insensitive to the p 
parametrization chosen; the value we find is —3500 ± 300 MeV^. This insensitivity is sig- 
nificant, for the p mass varies by some 10 MeV over the same set of the parametrizations. 
It is likely the consequence of the narrow uj width; in particular Fp/r^j ~ 20 [Q]. The s- 
dependence and phase of Ilf^{s), as per Eq. ( pi] ) and Eq. are also relatively insensitive 
to the p parametrization; this is shown in Tables |I| and |nT| . We conclude that U'p^ and 
Ul^iml) are 0.03 ± 0.04 and -300 ± 300 MeV^ respectively. The errors are such that both 
n' and n^, ,(m,^,) are consistent with zero. This is consistent with the conclusions of Costa 
de Beauregard et al. 0, and thus their implicit assumptions would seem to be justified in the 
resonance region. Note that the value of Re{Ilpa;(m2,)} continues to be —3500 ± 300 MeV^ 
in the presence of s-dependent or imaginary contributions to npi^(s); this is plausible as 
these corrections are themselves consistent with zero. Moreover, the inferred p parameters 
are insensitive to the inclusion of these effects as well. 

We also consider fits based on the form used by Barkov et al, Eq. (^2]). Fit E uses 
Particle Data Group values for the higher p resonances |jl[]; fit F uses the values adopted by 
Barkov et al. |T^. These parameters, along with the results of fits D, E, and F, are given in 
Table |V[ The values of np(^(m^) for fits E and F are determined from the fit to Eq. (p2D 



via the relation 

fr I'oTnZN _ ~3(m^ + dujmujTui)au) ,„_s 
I + a^j + ap' + Op" 

where follows from Eq. ([T3|) with mp,Tp m^,r^. The final error in the ITptj(m^) of 
fits E and F is determined by the full error matrix from the fit with Eq. (0). The form of 
Eq. ( p2|) is quite different from that of Eq. (|20|) , yet the values of Up^^^m'^) are comparable, 
cf. -3600 ± 300 MeV^ with -3500 ± 300 MeV^ from fits A-D. We favor the latter value as 
the p' and p" contributions in Eq. ( p2|) serve as sources of phase beyond that generated by 
Fp{s). We will examine whether such a construction is explicitly at odds with vr-vr scattering 
data. Note, however, that the scattering lengths extracted from fits D— F are reasonably 
close to the empirical value, though as s approaches threshold the impact of the imaginary 
part of the p' and p" contributions is minimized. 

The strength of p-uj mixing is commonly extracted from the ^ 2vr branching ra- 
tio [p|,p^, so that it is useful to reexamine this analysis and compare it with the current 
results. The physical \uj) and |p°) are related to the isospin-pure states \uji) and |p5) via the 
transformation 

(26) 

to leading order in isospin violation. Assuming that the vector mesons couple to conserved 
currents, ei and 62 are determined by requiring that the physical mixed propagator D^^i^s) 
has no poles, so that £1 = Ilp^{mlj — im^T^)/{m'^ — -|- i{mpVp — m^F^)) and 62 = 
^pu{ml - irripTp)/ (m^ - + i{mpTp - m^F^)). Then 
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/ + -I \ / + -I \ I np(^(m^ irripTp) -i o\ in'7\ 

TT + TT = TT + TT \uJl) + — FT^^ \Pll ^T) 

^""^""'^ (tt+tt-Ip?), (28) 



so that 



ml-m?p + i{mpVp - m^F, 



r(p ^ TT+vr ) plml \ {m? - )2 + {mpVp - m^T^y 



(29) 



where the overall p^^rn'^p/ (p^prri^) factor arises from treating {n^7i^\p'j), e.g., as the (7^^,^ cou- 
pling as per Ref. p9|. We assume that T{p — > vr+vr") = Tp and that (7r+7r~|p) = (7r+7r~|p/), 
as corrections to the latter are of nonleading order in isospin violation. Using the Particle 
Data Group values for the u —>■ 27t branching ratio, 2.21 ± .30%, and for the p and u reso- 



nance parameters |I| yields a central value of Iflp^;! = 3860 MeV , whereas using the Barkov 
et al. uo ^ 271 branching ratio, 2.3 ± 0.4 ± 0.2%, and rup = 775.9 MeV and Tp = 150.5 
MeV as per their analysis [0 yields Iflp^^l = 3950 MeV^ [^. Note that their u ^ 2% 



branching ratio is itself extracted from Eq. (0), though this and the PDG value use the 
estimated leptonic widths r(p — > e'^e~) and r(u; — e"^e~) to infer the ratio of the to p 
electromagnetic couplings, rather than adopt the SU{3) value of 1/3 as we do |31|]. This 
explains why the above results are of slightly larger magnitude than the result for Ilp^, 
-3500 ± 300 MeV^ which emerges from the direct analysis of e~^e — >■ tt+tt data in the 
time-like region. We prefer the latter analysis, however, for Ilp^ is accessed directly and is 
not subject to uncertainties in the p parametrization and its associated parameters. 

As the value of the p mass extracted by the analysis of Barkov et al. has been called 
into question on statistical grounds P^JI[], it is prudent to examine the sensitivity of our 
results to the chosen data set. Previously we have used the 82 data points with ^/s'^923 
MeV compiled by Barkov et al. and used in their analysis [l^. We have repeated fits A— C 



on the 40 points of the 1978 world data |T^ used by HL |]T5[ and on the 61 data points from 



the OLYA and CMD detectors reported by Barkov et al. [|T2[ in the s region of interest. 
Our results are shown in Table |V[ Note that the fit B errors refer to the errors generated 
by uncertainties in the data and by the error in the scattering length, respectively. We see 
that the parameters obtained for a given fit are consistent within errors for the three data 
sets. Note, moreover, that the striking insensitivity of np^(m^) to the p parametrization is 
manifest in all the data sets considered. 

The fits we have considered support a wide range of p masses and widths, rup, for 
example, ranges from 763.1 to 774.2 MeV, whereas Fp ranges from 144.0 to 157.0 MeV. The 
median values are 768.7 and 150.7 MeV respectively, remarkably close to the values reported 
by the Particle Data Group im . It is worth noting that the two-parameter Gounaris-Sakurai 



form, realized as either fit C, D, E, or F, recalling Eqs. ( pO| , p2D and Tables | and |IV|), 
consistently returns a value of 774 MeV, though fits D, E, and F all possess p widths some 
10 MeV smaller than that of fit C. Fits C and D only really differ in the manner they 
parametrize the s-dependence of the phase of the pion form factor, cf. h{s), Eqs. (p!0| ,|9|) , of 
fit C with g{s), Eqs. (|TB|,[T3), of fit D. 

The phase constraint can still be used to discriminate between the fits, for the phase of 
the form factor must be numerically that of / = 1, J = 1 scattering, if the latter is elastic. 
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We explore this issue in Fig. 0, in which the phase of fits A, C, E, and F are shown as a 



function of the invariant mass q, q = y/s, and compared with data |21|. In plotting the / = 1, 
1 = 1 phase shifts, we have adopted the energy-independent analysis of Hyams et al. |21[] and 
the energy- dependent analysis of Protopescu et al. The above analyses must assume 
isospin symmetry to separate the 5/ phase shifts, so that we omit the uj contribution from 
our plots of the phase of Ft^{s) determined from fits to the time-like pion form factor data. 
We have omitted fit D as its phase is essentially identical to fit C, which is shown. Fits C 
and D differ in how they parametrize the s-dependence of the phase of the pion form factor, 
yet this has no impact on the extracted tt-tt phase shift. The p widths for fits C and D do 
differ by some 10 MeV, but this is of no consequence for the accompanying phase shifts, 
as tan~^(ImF^/ReF^) — > 7r/2 as s m^. Fits B— F are reasonably consistent with data, 
though to judge this in detail we have computed the to the vr-vr scattering data for the 
fits shown in Fig. ^ 

xVdof = 1500/34 [A], 140/34 [B], 140/34 [C], 590/34 [E], 230/34 [F] (30) 

With such x^'s it is unlikely that any of the time-like pion form factor fits are truly correct 
models of the vr-vr phase shift data, though the employed data sets are themselves not 
always consistent within la. Note, moreover, that points in the p-uj interference region have 
been included in the computation, though in this region the u could influence the phase 
extracted from the p component of the form factor. Fit B is decidedly better than fit A; 
note that B differs from A only in that it incorporates the empirical scattering length. Fits 
C, E, and F are also similar in structure; they differ in the manner they parametrize the 
s-dependence of the form factor phase — fit C uses /i(s), Eqs. (]IO|,P), whereas fits E and F 
use g{s), Eqs. (|ll,0) — and in that the p' and p" contributions of fits E and F, Eq. (p2D and 



Table ^\ admit sources of phase beyond the p-uj sector. Fit C describes the vr-vr phase shift 
slightly better than fits E and F. As fits C and D yield phase shifts of comparable shape 
— the x^/dof for fit D is also 140/34 — it is the presence of the p' and p" contributions in 
fits E and F which is likely responsible for their poorer agreement with data. In particular, 
their inclusion increases the phase shift at large q. We have assumed that 5\ is purely elastic 
in the regime shown, though the measured elasticity vi\ does differ slightly from unity above 



900 MeV |]2l|. The structure of fits E and F suggests that this assumption be examined. If 
the phase shift is not strictly elastic, that is, if rj is not exactly unity, then the phase of the 
form factor is related to the scattering phase shift 5 via 

ImF (1 — r7cos2(5) ,„ . 

tan</) = = ^ . (31) 

^ ReF r/ sin 25 ^ ' 

If (1 — ■?]) <^ 1, then (j) and 5 are related by 

(l_^ (32) 

2,7] sm 



to leading order in (1 -r]). As r] < 1 and 5l> 100° for q > 800 MeV Eq. (H) implies 



that (j) < 6 as well. Thus the structure of fits E and F would seem inconsistent with the 
manner in which they fit the phase shifts, for the phase of fits E and F exceeds that of C 
for 5^900 MeV, at odds with the constraint of Eq. (|32|) . We conclude that the comparison 
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with the measured phase shifts indicates that fits B,C, and D would seem to be preferred. 
The average p parameters which emerge from these selected fits are mp = 773 MeV and 
Tp = 153 MeV. 

There is one last form factor constraint we can consider, for the pion form factor is known 



as s ^ from chiral perturbation theory. Through two-loop order [3S 



1 



(33) 



where {r )y is the electromagnetic charge radius of the pion squared and Cy is a low-energy 
constant. Note that fyij^) is the genuine loop contribution. A three-times-subtracted 
dispersion integral relates it to / = 1, / = 1 vr-vr scattering at the tree and one- loop level 
and to the pion form factor at the one-loop level ||38|; a convenient analytic representation 



of fvi-T^) is given in Ref. fy{-^) possesses contributions at both 0(s) and O(s^), 

yet the separation of the polynomial and dispersive pieces given in Eq. ( ^3]) suggests that 
its contribution to is very small — indeed, this is numerically the case |^8[. Crossing 

symmetry dictates the form of the 0{s) term and thus offers a consistency check of the space- 
like and time- like region data, (r^)^ has been measured in the space-like region by Amendolia 
et ai; adopting a fit in which \F^{t)\'^ = n/{l —t ■ |(r^)y)^ and -^,^(0) 



n 



1.000 ±0.009 

yields (t'^)^ = 0.431 ± 0.010 fm^ The data of Amendolia et al. has been reanalyzed by 
Colangelo et al. using the form dictated by chiral perturbation theory, Eq. (| 
find 



0.431 ± 0.020 ± 0.016 fm' 
3.2 ±0.5 ±0.9 GeV"^ 



and they 

(34) 
(35) 



where the errors refer to statistical and theoretical uncertainties, respectively. The agreement 
of with the value extracted by Amendolia et al. also indicates that the contribution of 
fy{-^) at 0{s) is small. The larger statistical error is a consequence of the two-parameter 
fit. 

In the time-like region, (r^)y is determined by fits to data above the two-pion threshold, 
s > Att"^. Ignoring the negligible Ilp^ contribution, one has for fits A and B, e.g., that 



{rTv = 6 A 




(36) 



whereas the (r^) contribution from F^^{s), Eqs. (p!l|,p!2 

Q / IT r ^ / ^2 



IS 



m?p ± dnipTp 




27rA;2 \^ ml) 



'nip ± 2kp 
2m^ 



1 r,m^- 



(37) 



Using the parameters of Table ^ and one finds 

{r% (fm^) = 0.30 ± 0.04 [A], 0.35 ± 0.03 ± 0.02 [B], 0.36 ± 0.03 [C&D], 0.40 ± 0.01 [E&F]. 

(38) 

Fit A does poorly on the a\ scattering length, so that its disagreement with the value of 
(r^)y from space-like data is not surprising. Fit B, however, is constrained to reproduce 
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a\ and is much closer to the Colangelo et al. value [^]. Note that the errors associated 



with fit B are to be added in quadrature; the first error arises from the error in the fits, the 
second from the error in a\. The agreement of fits C and D is comparable to that of fit B, 
whereas the results of fits E and F are within error of the Colangelo et al. value, when its 
error is taken into account as well. Thus, we see no real evidence for disagreement with the 
space-like data, and its consideration does not serve to distinguish the fits B— D, preferred 
by the phase shifts. 

A plurality of conclusions exists in the literature concerning the consistency of the space- 
like and time-like data Of those who have argued that they are inconsis- 



tent [p!5| , p| , the difficulty seems to be with data at larger t p!5|P^ , well beyond the region 
from which the charge radius is extracted. 

Let us conclude this section by summarizing our results for the p mass and width. Fits 
A— D are all distinct in character and generate excellent fits to the time-like pion form 
factor data. Presuming that the central values are normally distributed, we can compute 
the variance of the central values to infer the theoretical systematic error in a particular 
parameter. Using the m prescription, so that the form factor is given by Eq. as s — > m^, 
averaging over the results of fits A— D yields, 

= 770.5 ±1.1 ±4.4 MeV , = 153.5 ± 0.9 ± 4.0 MeV , (39) 

where the errors reflect the experimental statistical and systematic uncertainties and the 
theoretical systematic uncertainties arising from the p parametrization chosen, respectively. 
Note that we have discarded the error arising from the use of the empirical scattering length 
in Fit B in the above, as realistically one would refit the parameters should the scattering 
length be varied. These values differ from those of Barkov et al., irtp = 775.9 ± 0.8 ± 0.8 
MeV and Tp = 150.5± 1.6±2.5 MeV |T^, though the most significant difference is in our 
respective estimates of the model errors, cf. 4.4 MeV with 0.8 MeV for the p mass and 4.0 
MeV with 2.5 MeV for the p width. The theoretical systematic error we report is sensitive 
to the manner in which we determine the p mass and width. If, alternatively, we compute 
the mass and width from the value of the complex pole Sp, namely Sp = — imT, we find 
for fits A— D that Trip ranges from 756 to 757 MeV and Fp ranges from 141 to 142 MeV. A 
similar parametrization insensitivity of the mass and width defined from the complex pole 
position has been noted in Ref. 0. Note, however, for this latter choice, that the phase 
shift would not be 7r/2 at s = and that the resultant form factor near s = frfp would 
not be of Breit-Wigner form. Determining the mass and width as we have previously, but 
adopting the m prescription, noting Eq. (^, fits A— D yield 

mp = 774.3 ± 1.1 ±4.4 MeV , f p = 152.7 ± 0.9 ± 3.9 MeV . (40) 

If we discard fit A for its poor fit to the vr-vr phase shifts and repeat the above procedure 
for fits B— D only, we find 



and 



773.0 ± 0.7 ± 1.2 MeV , Fp = 153.3 ± 1.1 ± 4.6 MeV (41) 



mp = 776.8 ±0.7 ± 1.1 MeV , fp = 152.6 ± 1.1 ± 4.5 MeV . (42) 
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This step-wise elimination procedure may not be reasonable, however, as slight adjustments 
of fit A may yield an acceptable simultaneous description of the / = 1, J = 1 vr-vr phase shift 
and time-like pion form factor data with comparable values of the p resonance parameters. 
We favor our determination based on the time-like form factor data alone, Eq. (|39|) . 



IV. CONCLUSIONS 

We have determined the magnitude, phase, and s-dependence of the effective p-u mixing 
matrix element 'n.p,_j{s) from fits to e^e~ tt^tt~ data in the context of a theoretical 
framework which satisfies the analyticity, normalization, and phase constraints. We have 
considered a variety of descriptions which satisfy these theoretical constraints and have 
found the p-uj mixing parameters of interest to be insensitive to the manner in which the 
p is parametrized. Empirically <^ Tp, and this likely drives the above result. Recalling 
Eqs. (|2|J2|), we find 

Ilp^{ml) = -3500 ± 300 MeV^ 

^Ui^l) = -300 ± 300 MeV^ (43) 
n' = 0.30 ± 0.40 MeV^ , 

where n^^(m^) denotes the imaginary part of the effective mixing matrix element at s = 
and li'p^ characterizes the s-dependence of Iip^{s) about s = m^. Both the phase and s- 
dependence of np^(s) are statistically insignificant. It is not that these effects are numerically 
trivial, but rather that they are poorly constrained by current e~^e~ — > tt+tt" data. 

The value of ripi^(m^) we extract differs slightly from that extracted from the c<j ^ 27r 
branching ratio, Eq. (||), that is, \Ilp^{ml)\ ^ 3900 - 4000 MeV^. T{uj 2tt) is itself 
extracted from the time-like pion form factor data, and its connection to nptj(m^) explicitly 
involves the p parameters, which are relatively uncertain. It thus seems more appropriate 
to extract Up^^rn^) directly from e~^e~ data, as we have done. Were we to use Eq. p^ ) to 
extract B(a; — *• 27r) from our value of Ilp^(m^), not only would it be smaller than that which 
Barkov et al. reports [|T2|, but it would also be explicitly sensitive to the p parametrization. 



We have also systematically explored the p parameters associated with our parametriza- 
tions of the pion form factor. We find that the time-like data supports a range of p pa- 
rameters. Adopting the prescription currently favored by the Particle Data Group we 
find 

mp = 770.5 ± 1.1 ±4.4 MeV , = 153.5 ± 0.9 ± 4.0 MeV , (44) 

where the errors refers to the empirical and theoretical systematic uncertainties, respectively. 
It is worth noting that the differing prescriptions for the mass and width plague the direct 
comparison of the results of different groups, for the most recent world average for the 
p mass [|I| mixes results defined under different conventions The convention favored 



by the Particle Data Group has itself changed with time, c/. the m convention, noting 
Eq. (^, of Ref. |^ with the m convention, noting Eq. (^, of Ref. with no apparent 



adjustment of the reported values. We agree with the p width extracted by Barkov et al, 
Tp = 150.5 ± 1.6 ± 2.5 MeV, but their value for the p mass, rup = 775.9 ± 0.8 ± 0.8 MeV, 
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seems slightly large |T2 



It is the latter value that the Particle Data Group excludes on 
Note that our theoretical systematic errors are much larger than 
n, cf. 4.4 MeV with 0.8 MeV for the p mass. We have 



statistical grounds |jl;14 
those reported by Barkov et al. 
also studied the p parameters which result from different choices of data sets and find that 
the parameters extracted from the data of Barkov et al. are merely more precise. We 
conclude, then, that the data of Barkov et al. [|l6l are not likely wrong, but rather that the 
theoretical systematic errors associated with the extraction of the p parameters are larger 
than previously estimated. 
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TABLES 

TABLE I. Components of Ft^{s) for each fit. Note that "— " means tlpuj{s) = np^(m^). 
Fit n{s) Fred fipUs) 



A 


Eq. (|T6D 


Eq. © 




A' 


Eq. ® 


Eq. ^) 


Eq. 


Ai 


Eq. (|16D 


Eq. dll) 


Eq. dH) 


B 


Eq. (|16D (c fixed) 


Eq. (0) 




C 


Eq. (P) 


Eq. dll) 




D 


Eq. (|ll|) 


Eq. dll) 





TABLE IL Results from fitting the data of the Barkov et al. compilation from threshold 
through 923 MeV with F^{s) as per Eq. (|0|) and Table |. Note that = 0.13957 GeV. In fit B, 
the first set of errors associated with the p resonance parameters is determined by the statistical 
and experimental systematic errors in the time-like pion form factor data, whereas the second set 
of errors arises from the error in the input empirical scattering length. *Not a fitting parameter. 
Input. 



Parameter A B C A' A^ 



10^ • a (m-2) 


-2.80 ±0.88 


-0.672 ± 0.084 




-2.68 ±0.95 


-2.57 ± 1.06 


h 


0.24 ±0.52 


-1.008 ±0.027 


-1.403 ±0.023 


0.17 ±0.56 


0.10 ±0.63 


c {ml) 


11.4 ±8.1 


30.45* 


36.48 ± 0.66 


12.5 ±8.7 


13.7 ±9.6 


Pi (GeV-2) 


-0.65 ±0.25 


-0.27 ±0.12 


-0.16 ±0.12 


-0.69 ±0.25 


-0.57 ±0.27 


P2 (GeV-4) 


1.80 ±0.91 


0.76 ±0.46 


0.65 ±0.48 


1.85 ±0.88 


1.53 ±0.93 


P3 (GeV-6) 


-0.97 ±0.68 


-0.34 ±0.40 


-0.29 ±0.42 


-1.01 ±0.66 


-0.79 ±0.68 


fip^ (MeV2) 


-3460 ± 290 


-3460 ± 290 


-3460 ± 290 


-3500 ± 300 


-3480 ± 280 










0.027 ± 0.040 




Im fip^ (MeV2) 










-310 ± 280 


xVdof 


68/75 


68/76 


68/76 


67/74 


66/74 


Output 


Sp {ml) 


-12.46 


-125.7 




-14.21 


-8.214 


a\ (m-3) 


0.084 ± 0.043 


0.038 ± 0.002t 


0.0324 ± 0.0024 


0.079 ± 0.041 


0.073 ± 0.039 


rup (MeV) 


763.1 ± 3.9 


771.3 ± 1.3 ± 16 


773.9 ± 1.2 


763.7 ±4.1 


764.5 ± 4.5 


Tp (MeV) 


153.8 ± 1.2 


156.2 ±0.4 ±4.7 


157.0 ±0.4 


154.0 ± 1.2 


154.2 ± 1.3 


rfip (MeV) 


766.9 ± 4.0 


775.2 ±1.3 ±16 


777.8 ± 1.2 


767.5 ±4.2 


768.3 ± 4.5 


f p (MeV) 


153.0 ± 1.2 


155.4 ±0.4 ±4.7 


156.2 ±0.4 


153.2 ± 1.2 


153.4 ± 1.3 
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TABLE III. More results from fitting the data of the Barkov et al. compilation from threshold 
through 923 MeV with F^{s) as per Eq. (||) and Table |. Note that = 0.13957 GeV. Here the 
sensitivity of the s-dependence and phase of np^(s) to the p parametrization is examined. *Not a 



fitting parameter. 


^ Input. 








Parameter 


B' 


Bi 


C' 




10^ • a (m-2) 
b 

c (m^) 
Pi (GeV-2) 
P2 (GeV-4) 
/?3 (GeV-6) 
Hp^ (MeV2) 

Im iip^ (MeV2) 


-0.690 ± 0.087 
-1.002 ± 0.028 
30.43 lb 0.10 
-0.35 ±0.15 
0.93 ± 0.49 
-0.45 ±0.42 
-3510 ± 300 
0.034 ± 0.039 


-0.712 ±0.090 
-0.993 ± 0.030 
30.40 ± 0.10 
-0.25 ±0.12 
0.68 ±0.46 
-0.28 ±0.40 
-3480 ± 290 

-340 ± 270 


-1.408 ± 0.023 
36.62 ± 0.68 
-0.25 ±0.16 
0.82 ±0.51 
-0.40 ±0.43 
-3510 ± 300 
0.035 ± 0.039 


-1.413 ± 0.024 
36.80 ± 0.71 
-.14 ±0.13 
0.56 ± 0.48 
-0.22 ±0.42 
-3480 ± 290 

-340 ± 260 


xVdof 


67/75 


66/75 


67/75 


66/75 


Output 


mp (MeV) 
Tp (MeV) 
fhp (MeV) 
f p (MeV) 


-122.72 
0.038 ± 0.002t 

771.3 ± 1.3 ± 15 
156.2 ±0.4 ±4.3 
775.2 ± 1.3 ± 15 

155.4 ±0.4 ±4.3 


-119.01 
0.038 ± 0.002t 
771.5 ±1.3 ±15 

156.3 ±0.4 ±4 

775.4 ± 1.3 ± 15 

155.5 ±0.4 ±4 


0.0323 ± 0.0024 
773.9 ± 1.2 
157.0 ±0.4 
777.8 ± 1.2 
156.2 ±0.4 


0.0321 ± 0.0025 
774.1 ± 1.2 

157.0 ± 0.4 

778.1 ± 1.2 
156.3 ±0.4 
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TABLE IV. Results from fitting the data of the Barkov et al. compilation from threshold 



through 923 MeV with Ft^{s) as per Eq. ( p2D and Eqs. (11,12,13). However, note that fit D, as per 
Table |, uses Eq. (||) with Eqs. (|l|,|ll) in place of Eq. (H). For Fit D /3i = -0.16 ± 0.12 GeV-^ 



(52 = 0.65 ± 0.48 GeV-^ and [3^ = -0.29 ± 0.42 GeV'^. The values of tlpoj{ml) for fits E and F 
are determined from the fit to Eq. 



via Eq. (p5|). *Not a fitting parameter. ^Input. 







H; 


r 


mp (MeV) 
Tp (MeV) 

Op' 

mp/ (MeV) 
Fp/ (MeV) 
mp// (MeV) 
Fp// (MeV) 


773.9 ± 1.2 
146.9 ± 3.4 


774.2 ± 1.2 
145.7 ± 2.2 
(1.20 ± 0.34) • 10-3 
1.01 ±0.18 
-1.40 ±0.33 
1465t 
310t 
1700t 
235t 


773.8 ± 1.1 
144.0 ± 2.3 
(1.50 ±0.28) • 10-3 
0.467 ± 0.096 
-0.70 ± 0.20 
1290t 
200t 
1590t 
260t 


xVdof 


68/76 


74/77 


71/77 


Bp^ (MeV^) 
mp (MeV) 
f p (MeV) 
a{ (m-3) 


-3460 ± 290 
777.3 ± 1.2 
146.2 ±3.3 
(3.240 ± 0.060) • 10-2 


-3610 ± 310* 
777.6 ± 1.2 
145.1 ± 2.2 
(3.208 ± 0.036) • 10-2 


-3580 ± 310* 
777.1 ± 1.1 
143.4 ± 2.2 
(3.182 ±0.038) • 10-2 



TABLE V. p parameters and Ilf^{m'^) resulting from fits A— C of Table | and Eq. ( [20| ) to the 
time-like pion form factor data, \FT,{q^)\'^ . "40" denotes the 40 points of the 1978 world data. "60" 
denotes the OLYA and CMD data of Barkov et al., whereas "82" denotes the data set compiled 
by Barkov et al. 



Fit 


■nip 


Tp 


Upo; 


xVdof 


A40 


768 ± 12 


155.3 ±3.5 


-2970 ± 690 


41/33 


A61 


759.7 ±4.1 


152.8 ± 1.2 


-3340 ± 300 


38/54 


A82 


763.1 ±3.9 


153.8 ± 1.2 


-3460 ± 290 


68/75 


B40 


769.7 ± 3.6 ± 15.5 


155.7 ± 1.1 ±4.5 


-2970 ± 690 


41/34 


B61 


771.7 ± 1.3 ± 15.6 


156.3 ±0.38 ±4.6 


-3310 ± 300 


43/55 


B82 


771.3 ± 1.3 ± 16.1 


156.2 ±0.37 ±4.7 


-3460 ± 290 


68/76 


C40 


772.4 ± 3.6 


156.5 ± 1.0 


-2970 ± 690 


41/34 


C61 


774.0 ± 1.2 


157.2 ±0.4 


-3310 ± 300 


43/55 


C82 


773.9 ± 1.2 


157.0 ±0.4 


-3460 ± 290 


68/76 
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FIG. 1. The absolute square of the time-hke pion form factor, li^vrCs)]^, plotted versus the 
invariant mass q of the 7r"'"7r~ pair. Fits A (solid line) and B (long-dashed line), as per Eq. (20) 
and Table |, are shown with the data compiled by Barkov et al. The A fit, as per Table |, of Heyn 
and Lang (dashed line) is shown for reference. 
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FIG. 2. The 1 = 1,1=1 tt-tt phase shift 5j extracted from the phase of the time-hke pion form 
factor as determined by fits A (sohd hne), B (dotted hne), C (dashed hne), E (long-dashed hne), 
and F (dot-dashed hne) , given by Eqs. (|o|, ||) and Tables |, 0, and 0, plotted versus the pion-pair 
invariant mass q, along with the data from Ref. [26]. Note that the p-uj mixing contribution to 
the time-like pion form factor phase has been omitted, to facilitate comparison with the empirical 
phase shifts. 
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